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David S. Dean/ Pierre Le Doussal/ Satya N. Majumdar/ and Gregory Schehr^ 

^ Univ. Bordeaux and CNRS, Laboratoire Ondes et Matiere d’Aquitaine (LOMA), UMR 5798, F-33400 Talence, France 
^ CNRS-Laboratoire de Physique Theorique de I’Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex, France 
^ Univ. Paris-Sud, CNRS, LPTMS, UMR 8626, Orsay F-91405, France 

The ground state properties of N spinless free fermions in a d-dimensional confining potential are 
studied. We find that any n-point correlation function has a simple determinantal structure that 
allows us to compute several properties exactly for large N. We show that the average density has 
a finite support with an edge, and near this edge the density exhibits a universal (valid for a wide 
class of potentials) scaling behavior for large N. The associated edge scaling function is computed 
exactly and generalizes to any d the edge electron gas result of Kohn and Mattsson in d = 3 [Phys. 

Rev. Lett. 81 , 3487 (1998)]. In addition, we calculate the kernel (that characterizes any n-point 
correlation function) for large N and show that, when appropriately scaled, it depends only on 
dimension d, but has otherwise universal scaling forms, at the edges. The edge kernel, for higher d, 
generalizes the Airy kernel in one dimension, well known from random matrix theory. 

PACS numbers: 05.30.Fk, 02.10.Yn,02.50.-r,05.40.-a 


Over the past few years, experimental developments 
in the construction of optical traps and cooling protocols 
have allowed the study of systems of confined utracold 
atoms mm- Non-interacting fermions or bosons at low 
temperatures are of particular interest as they exhibit 
purely quantum effects. For example, Bose-Einstein con¬ 
densation mm has been observed experimentally in sev¬ 
eral cold atom systems. There are also nontrivial quan¬ 
tum effects in non-interacting fermionic atoms, arising 
purely from the Pauli exclusion principle. A well studied 
example is a system of N spinless fermions in a one¬ 
dimensional harmonic potential, V{x) = 

At zero temperature, T = 0, the many-body ground 
state wavefunction - ,xn) can be easily com¬ 

puted from the Slater determinant yielding nnnun] 

|4'o(a:i, • • • , xn )? = — (1) 

ZN // 

where a = y^mojjh has the dimension of inverse length 
and Ziv is a normalization constant. The squared wave- 
function |T'o(a;i,--- ,Xn)\'^, characterizing the quantum 
fluctuations at T = 0, can then be interpreted as the 
joint distribution of the eigenvalues of an A^ x A^ ran¬ 
dom matrix belonging to the Gaussian Unitary Ensemble 
(GUE) [nilll]. Consequently, several zero temperature 
properties in this one-dimensional fermionic system have 
been computed analytically iMa using established re¬ 
sults from random matrix theory (RMT). For example, 
the average density of fermions is given, for large N, by 
the celebrated Wigner semi-circular law [13 E]: 

with sharp edges at ±'/2N/a. These sharp edges get 
smeared for finite, but large, N over a width wn ~ 
and the density near the edge (say, the right one) is de- 
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FIG. 1. (Color online) (a) Plot of pN{x,y) vs. x and y for 
N — 28 fermions in a harmonic potential in d = 2 obtained 
from the exact evaluation [T^ of Eqs. (j^ and ( |10[ ). (b) Plot 
of pN{x,y) in the {x,y) plane predicted from the asymptotic 
formula Eq. (|^ with U(x) = and p « fica[r(d -|- 

1) (c) Comparison of (a) and (b) for y = 0 (the region 

below the asymptotic result is shown via the shaded region). 


scribed by a finite size scaling form [Ml Hi! 
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with Wn = /{ay^) and Fi{z) = [Ai'{z)Y — 

z[Ai{z)]'^, where Ai(z) is the Airy function. Far to the 
left of the right edge, using Fi{z) ~_y/z|/7r as z ^ —oo, 
one can show that the scaling form (1^ smoothly matches 
with the semi-circular density in the bulk ([^. The edge 
scaling function F’i(z) has been shown [3 to be univer¬ 
sal, i.e., independent of the precise shape of the confin¬ 
ing potential V{x). In addition, all n-point correlation 
functions at zero temperature can be expressed as deter¬ 
minants constructed from a fundamental quantity called 
the kernel (see later for its precise definition). For large 
N, and away from the edge, the appropriately scaled ker- 
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nel converges [3] to the universal sine-kernel form known 
from RMT [131 HI]- In contrast, near the edges, it ap¬ 
proaches [9] the universal Airy kernel, also well known 
from RMT [T71[T5]. Recently, for one-dimensional traps, 
several of the zero temperature results for the density as 
well as the kernel have been extended to finite tempera¬ 
ture [Ti^ . 

In many experimental setups, the optical traps are ac¬ 
tually in higher dimensions d > 1. For d > 1, there is 
no obvious relation between free fermions in a trap and 
RMT. Consequently, calculating analytically the zero 
temperature properties of spinless free fermions in a d- 
dimensional confining potential is a challenging problem. 
The bulk properties can be estimated rather accurately 
using the local density approximation (LDA, also known 
as the Thomas-Fermi approximation) For ex¬ 

ample, the bulk density for arbitrary confining potential 
V (x) is given, within LDA, by 


1 / m xd/2 - R(x)]‘^/2 

JV\2^J r(d/2 + l) 


( 4 ) 


where /r is the Fermi energy. This generalizes the 
semi-circular law ([^ in higher dimensions. Within 
LDA, one can even estimate accurately the two-point 
density-density correlation function within the bulk [see 
Eq. (18)]. However, near the edges (i.e., near the classical 
turning points where V (x) = fx ), LDA breaks down and 
the problem becomes more complex. Several studies have 
estimated the density and the correlation functions near 
the edge based on the extension of the semi-classical ex¬ 
pansion beyond the classically allowed regime [2n[23H26] . 
Notably, Kohn and Mattsson introduced an approximate 
method to calculate the density near the edge in d = 3 
[25| - they called it the edge electron gas. They provide 
a closed expression for the edge density in d = 3, but not 
for the edge correlation functions. For d ^ 1,3, neither 
the edge density nor the edge correlation function has 
been computed. In addition, unlike in d = 1, a system¬ 
atic analytical approach, valid both in the bulk as well 
as at the edges, is presently missing. 

In this Letter, we show that, in any dimension d, the 
n-point correlation functions in the ground state have 
a determinantal structure: it can be expressed as an 
nxn determinant det(x^,Xj) with 1 < i, j < n where 
Ar^(x, y) is the kernel, which exhibits universal scaling 
properties for large N (i.e., valid for a generic class of 
trap potential). This allows us to treat the bulk as well as 
the edge within the same unified approach and provides a 
rigorous basis for the edge electron gas in all dimensions. 
For simplicity, focusing first on the harmonic potential 
y(x) = where r = jxj, let us summarize our 

main results. First we recover the average density (one- 
point correlation) of the fermions at zero temperature for 
large N which is indeed given by Eq. with y(x) = 
and fi « hu![r{d + 1) is the Fermi energy 

for large N. The limiting density has a finite support of 
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FIG. 2. (Color online) Plot of the scaling functions Fd{z) in 
Eq. 0 for d = 1,2,3 (top to bottom) for the density near 
the edge. The oscillatory structure of the scaling function 
becomes less pronounced as the dimension d increases. 


radius redge = (v^/o^) where a = 

^JmujJh (see Fig. [^for d = 2). Our main result is to show 
that, in a region of width wn = bdN~^ (where bd = 
[r(d -b 1)]”®^ /{ay/2)) around the edge, the kernel takes 
a scaling form: iF^(x, y) « ^/Cedge 
with the scaling function 



( 5 ) 


where a„ = a • redge/j’edge and = b • redge/?'edge 
are projections of a and b in the radial direction and 
Aii{z) = Ai(u)du. Furthermore the scaling function 
/Cedge(a, b) is universal, i.e. valid for any spherically sym¬ 
metric potential V (x) [37] . Substituting a = b in Eq. (0 , 
we also obtain the scaling form of the edge density 

Pedge(x) = > (6) 

with 

1 

Fd{z)= 4 d d / duuiAi{u + 2^/^z), (7) 

r(§ -b l)2Tr7r2 do 

a plot of which is shown in Fig. [^ford= 1,2,3. Ind = 1 
it reduces to the RMT result Fi{z) = [Ai'{z)Y — z[Ai{z)Y 
mentioned earlier [EKn]. For d = 3 this result gives a 
rigorous proof, from first principles, of the edge electron 
gas result obtained in [35] . 

We start with N spinless free fermions in a d- 
dimensional potential y(x). The single particle 
eigenfunctions d’k(x) satisfy the Schrddinger equation, 
idd’k(x) = ek'0k(x), where H = —h^/{2m)'S/‘^ -b y(x) 
is the Hamiltonian and the energy eigenvalues Ck are 
labelled by d quantum numbers denoted by k. At 
zero temperature, the ground state many-body wave- 
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function can be expressed as an N x N Slater deter¬ 
minant, 4'o(xi, • • • , xat) = (l/-\/]V!) det[V'k(xj)], con¬ 
structed from the N single particle wavefunctions with 
energy up to the Fermi level /i [55] . For a sufficiently con¬ 
fining potential, fj, generically increases with increasing N 
PHD]. For example, for a d-dimensional harmonic oscil¬ 
lator y(x) = V{r) = muP'r'^/2^ « Hu![T{d + 1) 

Using det(A) det(i?) = det(Ai?), the squared many-body 
wavefunction can then be expressed as a determinant 

|^'o(xi,-- • ,x7v)P = ^ det (8) 

where the kernel iF^(x, y) is given by 

K^{x, y) = ^ 0{fi - ek)V’k(x)V'k(y) • (9) 

k 

Here 6{x) is the Heaviside theta function. As in d = 1 
[see Eq. 0]’ the squared wavefunction in ([^ can be in¬ 
terpreted as the joint probability density of N points in a 
d-dimensional space. By integrating out the N—n coordi¬ 
nates of |'I'o(xi, • • • ,xjv)p in Eq. 0, one can show (see 
Supp. Mat. [5^) that the n-point correlation function 
can be expressed as an n x n determinant whose entries 
are given precisely by the kernel iF^(x, y) in Eq. 0. In 
particular, for n = 1, the density /ojv(x) is given by 

Pn{^) = . (10) 

Thus, the knowledge of the kernel provides a complete de¬ 
scription of the statistical properties of the ground state. 

To compute the kernel, we first recall a well known 
relation between the kernel and the propagator of the 
single particle quantum problem |30j . Taking derivative 
of Eq. 0 with respect to (w.r.t.) /r and performing a 
Laplace transform w.r.t. ^ of the resulting relation, one 
finds 

9A:^(x, y) = e"* + '0k(x)V'k(y) • (H) 

dn ^ 

The right hand side of 0 is simply the single particle 
propagator G(x,y;t) = (x|e“ft^|y) in imaginary time. 
Integrating the left hand side of Eq. 0 by parts and 
inverting the Laplace transform, using Bromwich inver¬ 
sion formula, gives (see |30l |3T| and Refs, therein) 

^ ^ ^exp G(x,y;t) , (12) 



where F denotes the Bromwich integration contour. 


While Eq. (121 is general and holds for arbitrary poten¬ 
tial U(x), calculating the propagator G(x,y;t) explicitly 
is hard for general V (x). Below, we first focus on the spe¬ 
cific case of a harmonic oscillator for which G is known 
explicitly [55] 

/ 2 \ ^/2 2 


where Q{x, y; t) = (x^ -|- y^) cosh (w t) — 2 x • y. General 
potentials V (x) will be considered later. 

To extract the large N behavior of the kernel from Eq. 
(121, we need to perform a small t expansion of the prop¬ 
agator G(x,y;t) since fj, ^ is also large. Similar 

short time expansions have been carried out previously, 
especially in the context of nuclear [m [55] as well as 
chemical [34] physics. In this paper, we require such an 
expansion in general d as detailed in the supplementary 
material |29j . This expansion is however different from 
the Wigner-Kirkwood semi-classical expansion in powers 

of h [55] . 

Global density. We first evaluate the global density 
Pjv(x) in Eq. (10) by putting x = y in Eq. (12). The 
dominant contribution to the Bromwich integral in Eq. 
(12) with X = y comes from the small t region. Ex¬ 
panding the propagator to leading order for small t, 
the integral can be done explicitly to give the result in 
Eq. (lij) with U(x) = and fj. « /ia;[r(d-|- 1) 


The normalization condition f dx pj^ (x) = 1 fixes the 
Eermi energy p « huj[T{d + 1) The density thus 

has a radially symmetric finite support that vanishes 


at the edge as ~ (redge — (see Eig. 0, where 

J'edge = \/2^/(ma;2) ~ Since N particles are 

packed within a volume of radius fedge ~ 
typical inter-particle distance Gyp can be estimated very 
simply: ~ implying Gyp ^ 

Edge density. We next investigate the density near 
Tedge for finite but large N. To derive the asymptotic edge 
behavior, we again start with the propagator G(x,x;t), 
but now we set r = |x| = fedge + where (j) is yet 

to be determined and bd = [r(l -|-d)]~“/(a-\/2). Ex¬ 
panding the propagator for small t and keeping terms up 
to order 0{t^), we find that (p = l/(6d) in order that 
the two leading terms scale in the same way for large N 
with z fixed (see Supp. Mat. [5^). Subsequently, eval¬ 
uating the kernel in Eq. (12) and using Eq. (10), upon 
identifying wn = bdN~^/'^ , the edge density satisfies 
the scaling form in Eq. 0 where the scaling function is 
given by 

dr 1 


Fd{z) = 


^ -tz+tG12 

2TTi 


(14) 


Using the integral representation of the Airy funct ion . 
Ai{z) = l/(27ri) Jp dr the integral in Eq. (14) 

reduces to the expression announced in Eq. 0 (see Supp. 
Mat. [5D]). This thus generalizes to arbitrary d the Id re¬ 
sult, Fi{z) = [Ai'{z)]‘^ — z[Ai{z)Y, obtained from RMT 
[I6l[ll]. The asymptotic behaviors of Fd{z) can be com¬ 
puted explicitly (see Supp. Mat. [29]) with the result 


Fd{z) 


(8^)- 

(47r)- 


. tf+3 


^3/2 


as z 


CX) 


r(d/2 + i)' 

One can show that when z 


as 


— 00 


(15) 

(16) 


— 00 , i.e., when r ^ redge, 


the asymptotic behavior in Eq. (16) matches smoothly 
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with the bulk density in Eq. Q with E(x) = We choose D = 2^/^ and use this in Eq. (19). Using 

and ^ « hu}\r{d + 1 ) . 

Bulk kernel. We next consider the large N scaling be¬ 


havior of the kernel K^(x,y) in Eq. (12) where the two 


points X and y are both far from the edge, while their rel¬ 
ative separation |x—y| is on the scale of the inter-particle 
distance ~ . We start from Eq. (12) and 

in Eq. (13) we rewrite (3(x, y;t) = (x — y)^ + (x^ -|- 
y^)(cosh(wt) — 1). Expanding the propagator for small 
t to leading order we obtain 


/ W /" dt 1 

\27rtLiy 7r 27rz tf-i-i 


(^-V(|>c|))t o3(x-y)^ 


K^{x,y) 

\ ^iiuj / Jr 

(17) 

where U(|x|) = V{r) = muj‘^r'^/2. Fortunately, this inte¬ 
gral can be done exactly (see Supp. Mat. [13] )• We find 
that the bulk kernel has the scaling form, Er^(x, y) « 
- y|/£), where £ = [iVp 7 v(x) 7 d]“^/‘^ is the 
typical separation in the bulk and 7 ^ = 7 r‘^/^[r(d/ 2 -|- 1 )]. 
The bulk scaling function is given explicitly by 


^bulk(a^) = 


'^^ 72 ( 23 ^) 


(18) 


(Tra:)"^/^ 

where is the standard Bessel function of the first 

kind. In d = 1, using Ji/ 2 {z) = yj2 /(ttz) sinz, our result 
in Eq. ( |18| ) again reduces to the standard sine-kernel in 
RMT |13| . The result in Eq. (18) is in full agreement 
with the heuristic derivation using the EDA [ID] (see also 
Supp. Mat. [H]). However, the LDA becomes invalid 
near the edge where the local density is rapidly varying. 
We will see below that our approach yields exact results 
even in this edge regime where the LDA fails. 

Edge kernel. Turning to the large N behavior of the 
kernel Ar^(x, y) near the edge, we set x = redge + wn bl 
and y = redge + wn b, following the scaling of the 
edge density in Eq. (j^. Here redge denotes any point 
on the boundary of the support of the global density 
with |redge| = redge = 1 / 2 / 1 /(mw^) As be¬ 
fore, the width wn = bd with bd = [r(d -I- 

1)]“^/(®‘^)/(q;-\/ 2)- Thus a and b are dimensionless vec¬ 
tors. We substitute these scaling variables x and y in 
Q(x,y;t) = (x—y)^-|-(x^-|-y^)(cosh(a;t) — 1 ) and expand 
Q up to order for small t. Substituting these results 


in Eq. (12) and after a suitable change of variables (see 
Supp. Mat. [23]), one arrives at the edge kernel 

'■ dr 


Kt,{yi,y) 


1 


1 


CdW% /r 27ri r2+i 




{o-n.+bn)T 

21/3 


(19) 

with Cd = ( 23 yDr)'^, and where an = a ■ redge/redge and 
bn = h ■ redge/redge are projections of a and b in the 
radial direction. One can make a further simplification 


of Eq. (19) by using the integral representation of the 


diffusive propagator 


(a-b)2 


( 47 rZ? t ) 2 


(27r) 




V —zq- (a—b) 


( 20 ) 


subsequently the integral representation of the Airy func¬ 
tion Ai(z) mentioned earlier, we arrive at the main for¬ 
mula given in ([^ for the scaling behavior of the edge 
kernel. Putting a = b in Eq. (§, followed by an in¬ 
tegration by parts, one can check that /Cedge(a, a) re¬ 
duces to Ed(|a|) in Eq. ([^. Also one can verify, after 
a few steps of algebra (see Supp. Mat. [23]), that in 
d = 1 Eq. (j^ reduces to /Cedge(aj&) = Kxhy{a^b) = 
{Ai{a) Ai'(b) — Ai'(a) Ai{b))/{a — b) is the standard Airy 
kernel [13 [^. 

General potential. One naturally wonders to what ex¬ 
tent these results for the harmonic oscillator are univer¬ 
sal, i.e., hold for more general potentials U(x). We first 


note that for general U(x), Eq. (12) still holds, though 


the Fermi energy /i and the propagator G(x, y; t) depend 
on U(x). The dependence of the Fermi energy p. on N 
can be easily estimated for large N using semi-classical 
approximation PHD]. In contrast, G(x, y;t) is hard to 
compute for general U(x). However, as discussed after 
Eq. (12), for large N, we only need the small t expansion 


of G in general dimension d for arbitrary potential V (x). 
Using the results from [23], we find that to leading order 
in t the global density is given by Eq. Evaluating 


Eq. (12 1 near the edge, as in the harmonic oscillator case. 


we find that the edge density is again given exactly by 
Eq. (§, where only the location of the edge redge and 
the width wj^ depend non-universally on U(x), but the 
scaling function Fd{z) is universal, i.e., is the same for all 
spherically symmetric potential U(x) = V{r) ^7]. The 
bulk kernel is also given as in Eq. (18), the only U(x)- 
dependence appears in the scale factor yd used in defining 
the pair of dimensionless vectors u and v. Similarly, we 
find the edge kernel, appropriately centered and scaled, 
is also universal (see Supp. Mat. [33]). 

Conclusion. In this Letter, we have studied the ground 
state properties of N spinless free fermions in a d- 
dimensional confining trap. We have shown that the 
n-point correlation functions have a determinantal struc¬ 
ture for all d, with large N scaling forms that are uni¬ 
versal, i.e., independent of the details of the trap poten¬ 
tial. Our results recover the bulk properties predicted 
by the heuristic LDA. However, near the edge where this 
approximation fails, our method predicts new universal 
exact results in all dimensions d. In d = 1, we recover the 
results from RMT and in d = 3, our results recover rigor¬ 
ously the edge density result of Kohn and Mattsson [33 . 
Furthermore we provide an explicit expression in all d for 
the n-point correlation functions at the edge. Our results 
can be extended to finite temperature [35] in all dimen¬ 
sions. Finally, it would be interesting to see whether one 
can measure these universal edge scaling functions in ex¬ 
periments on cold atoms. 
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Supplementary Material 

We give the principal details of the calculations described in the manuscript of the Letter. 


DETERMINANTAL STRUCTURE FOR THE n-POINT CORRELATION FUNCTION 


Consider a situation where the joint distribution of N points |'I'o(xi, • • • ,xjv)P can be expressed as a determinant 


l^o(xi,--- ,X7v)P = ^ det A:^(xi,xP 
iV ! 


( 21 ) 


where the kernel ir^(x, y) has the reproducing property (which is the case for the kernel (8) in the Letter, as it is 
easy to check) 


J Kf,{x,z)Kf,{z,y)dz = Ki,{x,y) . 


The n-point correlation function i?„(xi, • • • ,x„ ) is defined as 

m 


ddn (x^l ; 


A' f 

= TTt -U / dXn+1 ■ ■ ■ dXN |^'o(xi, • • • ,XAr)|" 

(A - n)! J 


( 22 ) 


(23) 


(A-n)! 

obtained by integrating over N — n coordinates and keeping n coordinates {xi, • • • ,x„} fixed. If the kernel satisfies 
the reproducing property in Eq. (221, then there is a general theorem [T] that states that can be expressed as an 
n X n determinant 


i?„(xi, • • • ,x„) = det Kf,(xi,Xj) . 

l< 2 ,j<n 


(24) 


DERIVATION OF THE EDGE DENSITY IN EQ. (13) IN THE MAIN TEXT 

The single particle propagator for the d-dimensional harmonic oscillator is given by [5] 

d/2 


G(x,y;t) 


A- 


\ 27r sinh (w t) 


exp 


((x^ -I- y^) cosh (w t) — 2 X • y) 


2 sinh (w t) 

where a = y/mujIH. For the density [see Eqs. (9) and (11) in the main text], we need G(x,x;t), which is given by 

v2 \ 


(25) 


G(x,x;t) = - 


a 


27t sinh 


(ujt)) 


exp 


— tanh ( ^ 


where r = jxj. The kernel K^{x,y) is given by Eq. (11) in the main text, which reads 




(26) 


(27) 


where T denotes the Bromwich integration contour To evaluate the density /5jv(x) = (l/A)A^(x, x) for large N 
(equivalently for large p.) using Eq. 


we need the small t expansion of G(x, x;t) in Eq. (26). For small t, up 
to two leading orders, tanh(a;t/2) « ojt/2 — (wt)^/24. To compute the scaling behavior of the edge density, we set 
r = ?'edge + -z bd N~'^ where Cedge = \/‘^d‘l (ww^), bd = [r(d-|-l)]“3s j{aV2) and the exponent cj) is yet to be determined. 
We substitute this in Eq. (26) and use the resulting expression for G(x, x;t) in Eq. (27). Keeping terms up to order 
we get 


Kf,{x,x; t) 


( o? f dt 1 dLt 




(28) 


Note that the order 0{t^) term inside the exponential comes from the expansion of the prefactor [sinh (w 
The leading order terms pd/h cancel inside the argument of the exponential, leaving only three terms. Using /i « 
fujj\r{d -f we see that the three leading terms inside the exponential scale respectively as, Ti ~ 
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T 2 ^ and Ta ^ In the large N limit, we keep Ti ^ -j. ^ ^ Qf order 0(1). Then the next two 

terms scale as T 2 ~ ^ j^-i<t>-i/{ 2 d) ^3 then have two possibilities: (i) choose (j) = l/( 2 d) that 

keeps T 2 ~ ^ 0 ( 1 ), but then T 3 ^ N^/d^'i diverges as t 00, (ii) choose 0 = l/( 6 (i) that keeps ^ 0 ( 1 ) 

and T 2 ^ ]\f-‘^/(^d) ^2 vanishes as iV —>■ 00. Clearly the second choice, </> = l/( 6 d), is the correct one to make the 
scaling consistent. With this choice </> = l/( 6 d), one can also check that terms of order higher than O(t^) vanish as 
N —>■ 00. Hence, rescaling t = /{ijj[T{d + 1)]^/*^^'^^) r, one arrives finally at 

= (29) 

W% \ WN J 

where = hiiN~^, with bd = [r(l + d)] /{av2) and the scaling function Fd{z) is given by 

ft(.) = (4.)-/^/^ (30) 

Finally the edge density pedge(x) = (l/7V)A'^(x, x) is then given by Eq. (5) in the main text. 


DERIVATION OF EQ. (6) IN THE MAIN TEXT 

Here we show that Eq. (30) is identical to the expression given in Eq. ( 6 ) in the main text. We first use the identity 

1 1 r°° 


Tdl2+1 Y{d/2 + 1) 


dx 


in Eq. (30) to obtain 


Fd{z) = 


dx 


dr 

2711 


3 —t(x+z)+t ^/12 


r(d 72 +1) (4^)'^/2 jo 

After rescaling r —>■ and using the integral representation of the Airy function 


we get 




(31) 


(32) 


(33) 


(34) 


With a further rescaling x —>■ 2 Fd{z) in Eq. (34) reduces to the form given in Eq. ( 6 ) in the main text which 

reads 


Fd[z) = 


1 


r(f+ I)2w7r5 


nOO 

/ du Ai{u + 2^^^ z) . 

Jo 


In d = 1, the integral can be performed exactly. We start with the identity [4] 

and differentiate it twice with respect to z. Using the Airy differential equation Ai”(z) = zAi(z), one obtains 


) 

du y/uAi{z + u) = I"{z) — zl{z) = 7 r 2 ^/^ 


Ai' 






(35) 


(36) 


(37) 


(38) 


It then follows from Eq. (35), upon setting d = 1, that 

Fi{z) = At'\z) - zAi^(z) , 

thus recovering the well known RMT result 0 0 . One obtains similar quadratic forms in Ai{z) and Ai'{z) with 
polynomial coefficients in z in any odd space dimension by repeated application of the Airy operator (9^ — z) on /(z). 
For instance in d = 3: 


Fo^z) = -——(2z^Ai(z)^ — Ai(z)Ai'(z) — 2zAi'(z)^) 

1277 


( 39 ) 
















(40) 


In d = 2 one can use Airy equation and find 


F2{z) = —{-Ai'{23z) - 2izAii{2^z)) 

237r 

where Aii{z) = dxAi{x) also appears in Eq. (20) and can be expressed in terms of hypergeometric functions. 

Asymptotic behaviors of Fii{z). We first consider the z —)■ +oo limit. In this limit the Airy function has the leading 
asymptotic behavior [1] 




Ai{z 

Substituting this asymptotic behavior in Eq. ( |35[ ), expanding for large z, one gets to leading order 

Fdiz) « (Stt)”^ z~^ exp I —- " I as z —>• oo . 


(41) 


(42) 


Eor the other side z —>■ — oo, it is more convenient to use the representation in Eq. (30). We set z = — |z| and scale 
t|z| = t. This makes the order term to be |z|“^t^/12 which can then be dropped for large |z|. The resulting 
Bromwich contour integral can be easily evaluated to give the leading asymptotic behavior 


Fd{z) 


(47r)- 


r(d/2 + l) 


z 2 as z —>■ —oo . 


(43) 


DERIVATION OF THE BULK KERNEL 


In order to analyze the bulk kernel, it turns out to be convenient to rewrite the propagator G(x, y; t) given in Eq. 

d/2 


(25) in a slightly different form 
G(x,y;t) = 


27r sinh (w t) 


exp 


2 sinh (w t) 


((x- y)^ + (x^ +y2)(cosh(a;t) - 1)) 


(44) 


Substituting this in Eq. (27) and expanding the propagator to leading order in small t leads to Eq. (16) in the main 
text that reads 


where E(|x|) = V{r) = muj'^r'^/2. We then use the following integral representation [3] 


'=-y)2 


dt 1 


= (2v^) 


/r 27 rf \a 


(45) 


(46) 


where J,j(x) is the standard Bessel function of the first kind of index We use this identity to evaluate the kernel in 
Eq. (45). The result gets simplified upon replacing /i — Vir) in terms of the density using 


1 / m \d/'z [jji- 

Piv(x) « — r(d/2 + i) ■ 


(47) 


r(d/2 + l) 

The result can be cast in a simple scaling form 

~'^”'^^buik(|x-y|/£) (48) 

where I = is the typical separation in the bulk and 7 d = 7r'^/^[r(c?/2 + 1)]. The bulk scaling function 

is given explicitly by 

Jd/2 ( 23 ^) 


^bulk(3^) — 


(7ra:)‘^G 


(49) 


as in Eq. (17) in the main text, with Abuik(O) = l/7d- In d = 1, using Ji/ 2 {z) = ^J2/ (ttz) sinz, one recovers the 
standard sine-kernel /Cbuik(a;) = in RMT [I]. 
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It is important to note that this kernel is the Fourier transform of a Fermi-step function, using the formula: 


/|k|<fc/ 


(27r)'^ 




f kf 


\27r|x| 


d/2 


Jd/2ikf\x\) , 


(50) 


which can be verified in a straightforward way. To recover our result one must thus choose the local Fermi momentum 
as kf = kf{x) = 2lt = 2(iVpjv(x)7d)^^'^. This value is exactly consistent with the one obtained by the counting of 
states for a uniform system of density pAr(x), by setting N PAr(x) = — fc//(47r)‘^/^r(l + d/2). Thus, to 

describe correlations on scale i in the bulk, one can approximate locally the system by free fermions without any 
external potential, but at a fixed density pjv(x), assumed to be slowly varying on that scale. This then provides a 
more rigorous derivation of the results, within the bulk, obtained from the heuristic local density approximation. 


DERIVATION OF THE EDGE KERNEL 


In order to analyze the edge kernel we start with the propagator G(x, y; t) given in Eq. (44). Near the edge, we set 


X = Tedge + wn a and y = redge + wn b. Here redge denotes any point on the boundary of the support of the global 
density with |redge| = ’’edge = \/2.\i./(miiP-) ^ '/2/a\r(d + ^ section II for the edge density, we 

choose Wf^ = &d with bd = [r(d + /{a\/2). Thus denotes the width of the edge regime. The 

vectors a and b are thus dimensionless. As in the case of the edge density in section [[above, we expand the propagator 
for small t up to order 0{t^). Note that, for the edge kernel, we need to expand the propagator up to terms of order 
0{t^) (as opposed to the bulk kernel where it is sufficient to keep terms up to order 0{t)). As in section [| in the 
scaling limit, the term of order 0{t^) drops out for large N^ while the terms of order 0{t) and scale in the same 
way. Proceeding further as in section [| it is then straightforward to arrive at Eq. (18) in the main text that reads 


1 f dr 1 - _ (°Ti+i>n)T , 


(51) 


with Cd = (23 y^)'^ and where Un = a ■ redge/^'edge and = b • redge/''’edge are projections of a and b in the radial 
direction. One can make a further simplification of Eq. (51) by using the integral representation of the diffusive 
propagator 


e ■iOT 

{‘hvD t) 


d^q 


{2tt)° 


D q^T — ic[- (a—b) 


We choose D = 2^/^ and use this in Eq. (51). This gives the scaling behavior of the edge kernel, 


~ ^ ^edge 


W 


N 


X ^edge y ^edge \ 


Wn 


Wn 


where the scaling function is given explicitly by 


A^edge(a, b) — 


d<^q 


g iq.(a-b) f d.T 1 ^_(2-^/3^2^2-^^^(a„+b^))T+T^/3 


Defining Aii{z) as 


we arrive at 


(27r)‘^ Jy Sttz r 

K.,„(a,b) = j 


(52) 


(53) 


(54) 


(55) 


(56) 


The function Aii{z) defined in Eq. (55) can be actually expressed in terms of the Airy function itself. Indeed, 
differentiating Eq. 

oo) 


( [55[ ) with respect to 2 and using the integral representation of the Airy function in Eq. (33), one 
finds that Ai'i{z) = —Ai{z). Reintegrating with respect to z and using that Aii{z 


0, we get 


poo 

Aii{z) = / Ai{u) 

J Z 


du . 


(57) 
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REDUCTION TO THE STANDARD AIRY KERNEL IN d = 1 


Putting d = 1 in Eq. (56) we get 

^edge(^;^) — 


dq 


,iq(a-b) 


/22/3q2+2-l/3(a + &) 


Ai{z) dz . 


(58) 


Making a shift z = + 2 + b) + u gives 

/Cedge(a, 6) = / f Ai(u + + b)) du . 

J-oo Jq 

We next use a non-trivial identity involving Airy functions [4] 

r p. e-*9 iv-v') yij(22/3^2 ^ 2-^/^v + v')) = 2-% Ai{v)Ai{v') . 

J-oo 27r 

Choosing v = a + 2~‘^/^u and v' = b + substituting this identity in Eq. (59) and rescaling u —>■ 2~‘^/^u gives 

POO 

^^edge(o,^)= / du Ai{au) Ai{bu) . 

Jo 


(59) 


(60) 


(61) 


Since Ai{z) satisfies the differential equation Ai"{z) — zAi{z) = 0 we replace Ai{z) by Ai"{z)/z in Eq. (611. Next we 
use the identity 


1 


1 


{u + a){u + b) b — a [u + a u + b 


1 


1 


and integrate by parts. This then reduces Eq. (61) to the standard Airy kernel form 


/Codge(a, b) = ArAiry(a, b) = {Ai{a) Ai'{b) - Ai'{a) Ai{b))/{a - b) . 


(62) 

(63) 


KERNEL FOR GENERAL POTENTIAL 


We start from Eq. (27) [Eq. (11) in the main text] for the kernel where G(x,y;t) = (xje '‘'^jy) with the single 
particle Hamiltonian H = —S^/(2m)V^ + E(x). Expanding in the energy eigenbasis one gets 

G(x,y;t)=^vI;^(x)vI/s(y)e-^*/\ (64) 

E 

where 4'£;(x) satisfies the Schrodinger equation, —+ V e = E'^ e- It is then easy to verify that G satisfies 
the Feynmac-Kac equation 


dtG= ^V^G-\VG . 
2m h 


(65) 


To evaluate the kernel in Eq. (27) for large N (and hence for large fj,), we need the short time expansion of the 


propagator G(x, y;t) for general E(x). The Feynmac-Kac equation (65) is a good starting point for a perturbative 
small t expansion. We have computed explicitly the terms up to order 0{t^) using the perturbation theory of the 
Feynman-Kac equation [7]. Omitting details [7], we just quote the main results here up to order 0{t^): 


/ m \i I'm ,2l 




where 


'S'i(x,y)= [ duV{x+ {y -x)u), 

Jo 

<S' 2 (x,y)= [ duu{l-u){V^V){x+{y-x)u), 

Jo 

S^{x,y) = / du dv [min('u, ?;) — mi;] (Vy)(x-f (y — x) m) • (VU)(x-|-(y — x) w) 

Jo Jo 

— y— [ dw M^(l — m)^(V'^U)(x-I-( y — x) m) , 

4m do 


( 66 ) 

(67) 

( 68 ) 


( 69 ) 
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where (V^y)(x + (y — x) u) denotes the Laplacian evaluated at (x + (y — x) u) and similarly for (Vy)(x + (y — x) u) 
and (V‘*y)(x + (y — x) u). While the expressions for and S 2 have appeared in the literature before obtained by 
analyzing a discretized version of the path integral [ 8 ], the expression for S 3 has not been computed before, to the 
best of our knowledge. As a simple illustration, one can verify that for the Id harmonic oscillator, V(x) = 
one gets (setting for simplicity h = cj = m = 1 ): 


'S'i(x,y) = ^(x^ + xy) 

S 2 {x,y) = i 

S3{x,y) = + V + 7xj/) , 


which agree with the results from a direct expansion of the exact propagator given in Eq. (25) in d = 1. 


(70) 

(71) 

(72) 


Bulk kernel and global density. To compute the bulk kernel, we keep terms up to order 0{t) in Eq. ( 66 ). We 
substitute this expansion in Eq. (|27|) and use the integral representation in Eq. (46). This gives 


K^{x,y) = 


■m[y- Si{x,y)] 
27r^/i^(x — y)^ 


J d 


/ 2TO(x-y)2[^- 5'i(x,y7[ 


6»(/r-5'i(x,y)), 


(73) 


where J,^{z) denotes a Bessel function of the first kind with index i/. The support of the bulk part of the kernel is thus 
finite. The normalized particle density is given by pat (x) = (x, x) /N, and the Fermi energy fi is determined from the 

normalization condition f dxpjY(x) = 1 where the integral is over the support of the density. Using 5'i(x, x) = U(x) 
[see Eq. (67)] gives the global density for arbitrary potential 


1 / m \‘^/2 [p—U(x)]‘^/^ 

+1) ■ 


(74) 

[Npn(x)]-^/‘^, Eq. 


r(d/2 + l) 

At two generic points x and y in the bulk (far from the edges) with their separation |x — y| 

(73) simplifies to the scaling form 

Ar^(x,y) «-f"'^/Cbuik(|x-y|/.^) (75) 

where i = [NpN{x)yci]~^/'^ is the typical separation in the bulk and 7 ^ = 7 r^/^[r(d /2 + 1)]. The bulk scaling function 
is given explicitly by 


Abulk(^) — 


dd/ 2 {‘^x) 

(ttcc)'^/^ 


(76) 


as in Eq. (17) in the main text for the harmonic oscillator. The dependence on the potential V (x) enters only through 
the local density pn{x) and hence through the scale factor £. However the scaling function associated with the bulk 
kernel in Eq. (76) is completely universal for all U(x). 


Edge kernel. We now focus on the behavior of the kernel near the edge of the global density where the edge is defined 
by U(redge) = P- We consider (redge + Tedge + ^>0 where a' and b' are vectors such that ja'j, jb'j ^ [fedgel = ?'edge 
in a way we make more precise later. Henceforth, we focus only on spherically symmetric potential with a single 
minimum, e.g., U(x) = U(r) ^ r^, with p > 0, in the region r ^ Tedge (no other sc ale in the potential). For such a 


' r'P in the global density in Eq. (74) and using the normalization 
7 V 2 p/(d(p+ 2 )) ^g^j.gg jq_ Consequently, redge ~ S‘^£p+'^)) for 


potential, redge ~ Substituting U(x) = Vir) 

f dxpjY(x) = 1 , provides us with the estimate p ^ 
large N. 

To proceed, we set x = redge + a' and y = redge + b' in Eqs. ( 67|(69 ) and expand S'i,iS '2 and S 3 , assuming 
ja'I, jb^l ^ Tedge- As in the case of the harmonic oscillator, only the terms of order 0(t) and 0{t^) matter in the 
appropriate scaling regime. Keeping only these terms, the edge kernel is given by 


Ar^(x,y) 


dt 1 


\2'iTh/ yp27r*if+i 




exp 


m 

2 ht 


a! - b') 


exp 




(77) 


where VU = VU(redge) is evaluated at the edge. Here a’^ and denote the component of each vector normal to the 
edge. Following the analysis done for the harmonic oscillator, we then introduce the scaled dimensionless vectors a 
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and b defined via a' = wn a and b' = wn b, where the width wn has the dimension of length. Its precise value is 
fixed as follows. In order that both the terms of order 0{t) and 0{t^) in Eq. scale in the same way, we choose 
Wn as 


Wn = 


|W|2^ 

(2m)s 


(78) 


Note that for the potentials introduced above such that V{r) ^ rP (with p > 0), then |VE| ~ ?'edge implying 
Wn for large N. Thus for p > 1, the width shrinks with increasing N, while it increases for 

0 < p < 1. However, for any p > 0, wn ^ ^edge (recalling that redge A^ 2 p/(d(p+ 2 ))^^ Ppj. |.pg d-dimensional 
harmonic oscillator, p = 2, Eq. (78) gives wn = with bd = [r(d -|- I)]“^/(a-\/2) - in complete agreement 

with our previous analysis of the harmonic oscillator potential. Going back to the integral (0 we find that it is 
controlled by a time scale, which we denote by tAr, such that terms t an d are all of order unity, which imply 
respectively that tjv ^ ^ fully consistent with (78). We can thus estimate the magnitude 


of the neglected terms. The term 0(t^), from Eq. (68), scales as ~ it is thus small iff |vyp 


equivalently (p^^pp^ 1. If the potential does not contain other scales, this term scales as ^ ’’edge i® 

indeed negligible. Similarly it is easy to check that the neglected second term in S 3 in Eq. (69) is indeed small since 




..-p-2 


m|WP edge 

Substituting the scaling variables in Eq. ([77]) and rescaling time appropriately (following the same procedure as in 
the harmonic oscillator case in section 11 ) we get 


Ki,{x,y) 


CdW^ 


d Jr 


dr 1 


27ri • 


+1 


• e 28/3-r 


{o-n + bnh 

2 I /3 




(79) 


where Cj = (23y^)^. For the d-dimensional harmonic oscillator, one can verify that Eq. (79) reduces exactly to the 
result in Eq. (51). Repeating the analysis done after Eq. (51), we finally get the scaling behavior of the edge kernel. 


Kf,{x,y) 


d '^edge 


^ I’edge y ^^edge \ 


WN 


WN 




where the scaling function is given explicitly by 


/Cedge(a,b)= I (^2h^ + 


21/3 


(80) 


(81) 


where we recall that Aii{z) = Ai(u) du. Thus all the dependence on the potential E(|x|) is encoded in the width 
Wn of the edge regime, but the scaling function /Cedge(a,b) in Eq. (81) is universal, i.e., independent of 17(r). 


Edge density. Finally, putting a = b in Eq. (81) and using the definition pedge(x) = {l/N)K^{:x.,x) 
rc(^‘^/Cedge(a, a)/iV, we arrive at 


Pedge(x) « +a)) , 

N w% J (27r)‘' 


(82) 


(83) 


where a = |a| and Aii(z) = Ai{x) dx. To proceed, we rewrite Eq. (82) as 

pOO nOO 

Pedge(x)«— -— dqq^-^ Ai{z)dz, 

N {2TTWNr Jo d22/3(q2+o) 

where T^ = 27r‘^/^/r(d/2) is the surface area of the d—dimensional unit sphere. Making the change of variable 
u = < 7 ^ 2^/3 and integrating by parts gives us 




w 


N 


Wn 


(84) 


where the scaling function Fdiz) is given in Eq. ( [34| [which appears as Eq. (6) in the main text]. As in the case of 
the edge kernel, the dependence on H(|x|) appears only through the scale factor wn, but the scaling function for the 
edge density Fd(z) is universal. 
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